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The hydrodynamic description of transversally thermalized matter, pos- 
sibly formed at the early stages of ultra-relativistic heavy-ion collisions, is 
developed. The formahsm is based on the thermodynamically consistent 
approach with all thermodynamic variables referring to two-dimensional 
objects, the so-called transverse clusters, which are identified with the par- 
ticles having the same rapidity. The resulting hydrodynamic equations for 
a single cluster have the form of the two-dimensional hydrodynamic equa- 
tions of the perfect fluid. Since the clusters do not perform any work in 
the longitudinal direction, their energy is completely transformed and used 
to generate strong radial and elliptic flows that turn out to be compatible 
with the experimental data. 
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1. Introduction 

The data collected at RHIC indicates that matter produced in ultra- 
relativistic heavy-ion collisions behaves like an almost perfect fluid [HE]. 
This observation triggers many new developments of relativistic hydrody- 
namics of perfect and viscous fluids [HllllllEllTlElEllIOllIIlliailSllIllllll 
dH [iTl dSI ■ On the other hand, the hydrodynamic picture is still challenged 
by the two serious problems. The first one refers to an unexpected short 
thermalization scale that is required to describe the measured asymmetry of 
the transverse flow, the second one is connected with the difficulty to explain 
the measured HBT radii (the so-called HBT puzzle). 
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Recently, it was proposed that these difficulties can be avoided if at the 
early stages of the evolution of matter produced in ultra-relativistic heavy- 
ion collisions the hydrodynamic description applies only to the transverse 
degrees of freedom while the longitudinal expansion may be to large extent 
regarded as the motion of independent clusters [lOj. In this paper we explain 
in more detail and develop the concept of Ref. [19]. The idea of purely 
transverse equilibration has been analyzed previously by Heinz and Wong 
in Ref. [20] with the conclusion that it cannot be realistic since it does 
not lead to the large elliptic flow found in the corresponding 3-dimensional 
(3D) hydrodynamic calculations. Our conclusions differ substantially from 
those reached in [20] because of at least two reasons: Firstly, we compare 
the results of our model calculations to the present data rather than to 
other hydrodynamic calculations. Secondly, we use a different technical 
implementation of the concept of transverse thermalization and longitudinal 
free-streaming. 

The formalism introduced in [19] and developed in this paper is based on 
the thermodynamically consistent approach where all thermodynamic vari- 
ables refer to two-dimensional objects. We call them transverse clusters and 
identify with the particles having the same rapidity. The resulting hydrody- 
namic equations for a single cluster, i.e. for a fixed value of rapidity, have 
the form of the 2-dimensional (2D) hydrodynamic equations of perfect fluid. 
Since the clusters do not perform any work in the longitudinal direction, 
their energy is completely transformed and used to generate strong radial 
and elliptic flows that turn out to be compatible with the experimental data. 
In contrast to Ref. [20] our approach conserves entropy and, therefore, de- 
scribes perfect fluid. 

An attractive feature of our approach are short space and time scales 
characterizing the system at the moment when the substantial elliptic flow 
is formed. We expect that this feature may help to solve the HBT puzzle and 
the problem of early equilibration. The verification of this point requires, 
however, further developments, e.g. an implementation of the hadronization 
scheme which transforms partons into the observed hadrons. 

The paper is organized as follows: In Sects. [2]- [4] we introduce the phase- 
space distribution function, its moments, and the equations of transverse hy- 
drodynamics, respectively. The moments of the distribution function define 
the particle density current, the energy-momentum tensor, and the entropy 
current. The hydrodynamic equations are obtained from the conservation 
laws for the energy and momentum. In Sects. [5] and [6] we discuss our initial 
conditions and the Cooper-Frye formula. In Sects. [7] and 8 we present our 
results and conclusions. Finally, in the Appendix the technical details of our 
calculations are given. 
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Fig. 1. Visualization of the ansatz ([I]) - ([2]). Particles having the same longitudinal 
velocity form transverse clusters. Their dynamics is governed by the equations of 
2D hydrodynamics of perfect fluid. The rapidity distribution of clusters is deflned 
by the function /|| . 



2. Ansatz for the phase-space distribution function 



In our approach we assume the fonowing factorization of the phase-space 
distribution function f{x,p) into the longitudinal and transverse part, see 

Fig.m 

f{x,p) = f\\gcq{T,V,^±,^J_) , (1) 

where 

/|| = no<5(p||t - Ez) = no^^^^^. (2) 

Here the standard definitions of the energy and longitudinal momentum in 
terms of the rapidity have been used, 

E = m± coshy, = m± sinhy. (3) 

Similarly, the spacetime coordinates t and z have been expressed in terms 
of the proper time r and spacetime rapidity rj, 

t = T coshr/, z = T sinhry. (4) 

We note that the form of the distribution ([T]) - ([2]) is very much similar to 
that used by Heinz and Wong in Ref . [20] , the only formal difference is that 
the factor no/m± in Eq. ([T]) is replaced in Ref. [20] by the factor tq. With 
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our ansatz all thermodynamic variables which appear in the formalism are 
consistently defined as two-dimensional quantities and the resulting equa- 
tions describe perfect fluid. This is different from Ref. [20] where the authors 
treat their approach as non-equilibrium viscous dynamics which generates 
entropy. 

In Eq. ([2]) the dimensionless normalization parameter uq describes the 
density of clusters in rapidity. The role of this parameter may be compared 
with the role played by the time scale parameter tq introduced in Ref. [20j . 
Heinz and Wong interpret tq as the initial time for the hydrodynamic evo- 
lution. A closer inspection of their approach indicates, however, that this 
parameter plays a role of the global normalization having no relation with 
the initial time. The hydrodynamic equations of Ref. [.2QJ do not contain the 
To parameter. The hydrodynamic evolution may start at arbitrary value of 
the initial time, let us call it r = Tinit, and the final results are independent 
of the choice that has been made for Tinit- 

The equilibrium distribution g^q has the form of the two-dimensional 
equilibrium distribution function convoluted with the transverse flow. For 
simplicity we use the Boltzmann statistics and neglect the chemical poten- 
tial. In this case we have 

fifeq = exp , (5) 



where the transverse mass m_L is defined by the formula m± = ym'^ + P^- 
The transverse flow has the structure 

u'' = (uq, Ux, Uy,0) = {uq, 1?±,0), 

^0 - "^i = 1> (6) 

which explicitly underlines its transverse character (vanishing longitudinal 
component, Uz =0). 

3. Moments of the phase-space distribution function 

By calculating the appropriate momentum integrals of the distribution 
function one obtains the particle current A^'^, the energy-momentum tensor 
T'^'', and the entropy current S'^. Treating particles as massless with the 
ansatz ([T]) one finds 
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~ "° (27r)2 ^ ^ rmx ^'^'^ 



27rT 



(27r)2 rm_L 

= (9) 
2ttt 

Here we have introduced the degeneracy factor connected with the density 
of states in the transverse space. Treating our system as dominated by gluons 
we assume that Vg = 16. In Eqs. (I7| - ([9|) the 3D hydrodynamic flow 
has the structure 

= (coshr/MO) Ux, Uy, sinhr^uo). (10) 

We note that the correct normalization of the transverse flow u'^ leads di- 
rectly to the correct normalization of the four-vector U^, namely U^U^ = 1. 
We also note that the energy-momentum tensor includes an extra contri- 
bution which does not appear in the standard hydrodynamics. The extra 
term is proportional to the product V^V^, where the four vector has the 
structure 

V = (sinhr/, 0, 0, coshr/) . (11) 

The four vector is spacelike, V^V^ = — 1, and its origin is connected with 
special role of the longitudinal direction in our case, at r/ = it takes the 
form = (0, 0, 0, 1). With the presence of the term proportional to V^^V^ , 
the energy-momentum tensor is traceless, T'jj = 0, as required for massless 
particles. The matrix form of the energy-momentum tensor is given in the 
Appendix. 

The appearance of the extra terms in the energy-momentum tensor was 
already discussed in Ref. [2D]- We stress that our structure is simpler 
than that discussed by Heinz and Wong and this fact is connected with 
the identification of the two-dimensional thermodynamic properties of the 
system. With our ansatz for the distribution function, the time compo- 
nents of the currents (I7| - ([HI) reduce in the rest frame of the system to the 
two-dimensional thermodynamic densities divided by the proper time. In 
particular, in the frame where = (1,0,0,0) we find 

N' = ^n,, T^' = ^S2, S' = ^S2, (12) 
r r r 

where the appropriate two-dimensional densities are defined by the equations 

f d^P T^gT^ 
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S2 = -Vg J J^9eq ( In 5eg " 1) = (13) 

Those definitions may be supplemented by the definition of pressure, which 
is obtained from the thermodynamic relation 



e2 + P2 = Ts2. (14) 



A simple calculation gives 



The last equality, in agreement with the expectations for two-dimensional 
systems, yields the sound velocity 

cl = \. (16) 

We note that the factor 1/r in Eqs. (fT2]) is of pure kinematic origin and 
describes the decrease of three dimensional densities due to the increasing 
distance between the transverse clusters. At midrapidity we have dz = t dy 
and the rapidity densities per unit transverse area, dA = dxdy, are 

dN dE dS 

non2, -nrr = '^0^2, -tt-t = noS2. (17) 



dAdy dAdy dAdy 

Eqs. (flTl) are natural for 2D systems, and one can see that rapidity densities 
per unit transverse area change only if the temperature decreases, i.e., only 
if the transverse flow is present. We note that this type of behavior is 
completely different from the scenario assumed in the Bjorken model. 

The temperature dependence of our 2D thermodynamic variables is dif- 
ferent from the scaling of the energy density £3 ~ and pressure P3 ~ 
that was used in Ref. [Ml- The relations £3 ~ r^,P3 T^, and £3 = 2P3, 
when used in equilibrium, lead to the contradictory results for the value of 
the sound velocity. The first two yield S3 ~ and = 1/3, whereas the 
third one gives = 1/2. Consequently, the approach presented in Ref. [20] 
should be interpreted as an effective description of viscous dynamics, while 
in our approach the transverse dynamics is reduced to 2D hydrodynamics 
of perfect fluid with = 1/2. 

The large value of the sound velocity leads, as expected, to much stiffer 
than usual equation of state, and favors formation of the strong elliptic flow. 
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Fig. 2. Geometrical definition of the dynamical angle a. In our approach, instead 
of the two components of the transverse flow we use the magnitude of the flow v 
and the angle a. Both v and a are functions of r, r and 4>. 



The most important effect responsible for the formation of the strong flow is, 
however, the lack of the interaction between the transverse clusters. Since 
the energy of the clusters is not reduced by the work done in the longitudinal 
direction, it is exclusively transformed and used to generate strong radial and 
elliptic flows. 



4. Transverse hydrodynamics 

The hydrodynamic equations are obtained from the energy and momen- 
tum conservation laws 

^^T^'^ = 0, (18) 

where the energy-momentum tensor is given by Eq. jH]). We have checked 
that Eqs. (fTSl) imply the entropy conservation law 

^^,S^' = 0, (19) 

with the entropy current defined by Eq. ([H). The connection between Eqs. 
(fTSll and Eq. (fTOl) is the same as in the standard 3D hydrodynamics of 
perfect fluid where the entropy conservation law d^S^^ = is obtained by the 
projection of the energy-momentum conservation laws on the four-velocity 
of the fluid, i.e., from the formula UydyT^'^ = 0. 

The use of our form of the energy-momentum tensor in Eq. (fTSj) leads 
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to the following three equations 

d d d 

— {rs2UQ) + — {rs2Uovcosa) + — {s2UQVsma) = 0, 
OT or Ocp 

d d d 

— {rTuQv)+r cosa—{Tuo)+ sm a— {Tuq) = 0, 
OT or ocp 

^ f da V sin a \ dT cos a dT 

Tulv — + - sma— + — = 0. 

Var r J or r Oo 



(20) 



The first equation in (|20)1 is just the entropy conservation, compare Eq. (1191) . 
The second equation in ((20]) follows from the symmetric linear combination 
UidpT^^ + U2dfj_T'^'^ = 0, while the third equation results from the asym- 
metric linear combination U2d^,Tt'^ - C/i^^T^^ = q. In Eqs. ([20]) we used 
the cylindrical coordinates 



x2 + y2^ (/) = tan ^{y/x). (21) 



1 



The quantity v is the transverse flow, uq = {l — v'^) ^ , and a is the dynam- 
ical angle describing deviation of the direction of the flow from the radial 
direction, 

Vx = V cos{a + (p), Vy = V sm{a + (j)) , (22) 

see Fig. [2l The deferential operator d/dT is the total time derivative defined 
by the formula 

d d d vsina d 

— = ^ + vcosa—^ — . (23) 

dT OT or r od> 



The structure of Eqs. ()20]) is similar to the structure of the equations de- 
scribing the boost-invariant 3D hydrodynamic expansion, see [21] and [22], 

— [rs3Uo) + T:- {rs3Uovcosa) + — {s3Uovsma) = , 

OT or ocp T 

Q d d 

— {rTuQv) + r cos a— {Tuq) + sin a— (Tuo) = 0, 
OT Or Ocp 

^ o f da V sin a \ dT cos a dT 

TuqV — H sma— H — = 0. 

\dT r J or r dcp 

(24) 

They differ from Eqs. ([20]) by the presence of the term rs^uo/T on the 
right-hand-side of the first equation in ([24ll and by the different form of 
the entropy density. We recall that the three dimensional entropy density 
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S3 of massless particles with u internal degrees of freedom is given by the 
expression 

^3 = (25) 

The main physical difference between Eqs. ([20]) and (p4l) resides in the term 
{rs3Uo)/T leading to the decrease of the entropy and energy densities even 
in the case where the transverse flow is absent. The physical origin of this 
term is the presence of the work which is done in the longitudinal direction. 
In the case of the 2D expansion such a term is not present and the deposited 
collision energy is used only to produce the transverse flow. 

We note that ([20]) is valid even if the temperature T, the transverse flow 
V, its direction a, and the normalization factor no depend on the spacetime 
rapidity ry. The structure of the energy-momentum tensor ([8|) implies that 
all partial derivatives with respect to tj are multiplied by the expressions 
which vanish in the case rj = y. Consequently, Eqs. ([20]) are valid for any 
value of the rapidity, and they should be solved, with the corresponding 
initial condition, independently for each value of rj. This property shows 
explicitly that our system is not boost-invariant and may indeed be treated 
as a superposition of the independent transverse clusters. 

The change of the initial entropy density is nicely described by the global 
conservation laws which follow from Eqs. ([20]) and ([24l) . The integration of 
the first equation in ([20l) over the transverse spacetime coordinates yields 

oo 2lT 

J drr j d(j)T'^{T,r, (l))uQ{T,r, (p) = const. (26) 


On the other hand, for the 3D boost-invariant case one finds 

oo 2lT 

J drr J d<PT^ir,r,cl,)uo{T,r,^) = (27) 


For 2D hydrodynamic expansion, see Eq. ([26]) . the initial thermal energy 
may decrease only at the expense of increasing transverse flow. For 3D boost- 
invariant expansion, see Eq. ([27|, even without the transverse expansion the 
temperature drops, as is well known from the Bjorken model. 

We solve Eqs. ([20]) . and also for the comparison Eqs. ([24]) . using the 
technique presented in Ref. [22j which is a direct extension of the method 
proposed in Ref. [23]. This method satisfies very accurately the conservation 
laws ([26]) and ([27ll . We find this agreement as an important check of our 
numerical scheme. 
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5. Initial conditions 



The hydrodynamic equations (|20p are three equations for three unknown 
functions: T, f , and a. At the initial time r = Ti^a the transverse flow v 
is zero. We also set the angle a to be zero at r = Tinit- On the other hand 
the temperature profile at r = Tinit is not trivial and its asymmetry in the 
transverse plane generates the elliptic flow. 

Similarly to other hydrodynamic calculations we assume that the initial 
energy density at the transverse position point lf± is proportional to the 
wounded- nucleon density pwN at this point, namely 

£2 C^±) = oc PWN ( ■ (28) 

vr 

We note that the assumption (f28l) used by us for a 2D system is equivalent to 
the assumption S3 oc pwN used in 3D hydrodynamic codes. In our practical 
calculations Eq. ((28]) takes the form 



r(Tinit, 



PWN i'^. 



1/3 



(29) 



PWN (0) 

where the parameter Tj is the initial central temperature and the wounded 
nucleon density is obtained from the formula [24] 



PWN 



+ Ta 




a 
'A 



Ta 



A 



Ta 




+ X _L 



(30) 

In Eq. ([3Q]l a = 40 mb is the total nucleon- nucleon cross section and Ta (x, y) 
is the nucleus thickness function 



TA{x,y) = / dzp{x,y,z) . 



(31) 

Here p{r) is the nuclear density profile given by the Woods-Saxon function 
with the conventional choice of parameters used for the gold nucleus: 

-3 



Po 
a 



0.17fm" 
(1.12A^/3 -0.86^~^/^)fm, 
0.54 fm, A = 197. 



(32) 



The value of the impact parameter in Eq. (I30|l depends on the centrality 
class considered in the calculations. 
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Fig. 3. Geometric interpretation of the parameters used to define the shape of the 
hypersurface corresponding to a constant value of the final temperature, T = Tf. 



6. Cooper-Frye formula 

To calculate the transverse-momentum spectra at a certain value of the 
final temperature Tf we use the standard Cooper-Frye prescription 

fd^^p,'-^^g.,, (33) 



dy(Pp± (2vr)2 J Tm± 

where the hypersurface S is defined by the condition Tf = const. For 
cylindrically asymmetric collisions and midrapidity, y = 0, the transverse- 
momentum spectrum has the following expansion in the azimuthal angle of 
the emitted particles 

{l + 2v2{p±)cos{24>p) + ...). (34) 



dyd?p± dy 27rp± dpj 

Eq. ([34l) defines the elliptic flow coefficient U2, which may be calculated from 
([34l) as the asymmetry of the momentum spectrum, 

, ^ _ 1 fN{P±, (Ap = 0) - fN{P±, 4>p = f ) 

'''^^^^ " 2fM{p±,^p = 0) + fN{p±,cPp = ^)' 

(35) 

with /tv being a shorthand notation for dN/{dyd'^p±). 

We are of the opinion that for the essentially 2D expansion, the freeze- 
out criterion cannot involve the 3D energy density as it was proposed in 
|20| . Application of a 3D freeze-out criterion to a superposition of 2D sys- 
tems implies that the freeze-out is triggered, in the artificial way, by the 



^ Since we concentrate the evolution of partons this temperature cannot be directly 
interpreted as the freeze-out temperature. 
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increase of the relative distance between the 2D parts. With those remarks 
in mind we adopt a different strategy; we present our results for different 
final temperatures and check if the experimental data can be successfully 
reproduced for one of them. 

In our calculations we use the following parameterization of the hyper- 
surface S, 



t = d ((/), C, ry) sin C cosh 7/, z = c? ((/>, t?) sin(^ sinh??, 
X = d (0, C, ??) cos C cos y = d C, ry) cos C sin(/), 



(36) 



which also yields 



T = \/t'^ — z'^ = d ((/), C, ■?/) sin 

p = ^Jx^ + y^ = d(0,C,r?)cosC. (37) 

At any given value of the spacetime rapidity r/ the position of the point on 
the hypersurface is defined by the two angles, (j) and C, and the distance to 
the origin of the coordinate system, d{(j),C,v)^ see Fig. [3l The angle (j) is 
the standard azimuthal angle in the y — x plane, while the angle C is the 
azimuthal angle in the t — p plane. With the standard definition of the four- 
momentum in terms of the rapidity and transverse momentum, and with 
the standard definition of the element of the hypersurface dS^ in terms of 
the totally antisymmetric tensor namely 

dx° dx^ dx'^ , , , ,A /or^N 
dS, = e,„,, — — — d^d^dC, (38) 

we find the explicit form of the Cooper-Frye integration measure 



dT^aP^ = d^sinC 



d cos C {m,±_ sin Q cosh [r] — y) +p±_ cos C, cos {(p — (j)p)) 



dd 

-|-— cos C {—m±_ cos C cosh {rj — y) + p± sin ( cos {(j) — (pp)) 

dd . , , , , 
+7r7P± sm (0 - <pp) 



drjdcpdC. (39) 



dd ^ . , . . 
-|---— cot Q m_|_ smh (77 — y) 
or] 

Eq. ([37ll and the presence of the delta function in Eq. ([331) imply that 
the integration measure for massless particles appears always in a simple 
combination 
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rp± 



d cos C (sin C, + cos Q cos {(j) — (j)p)) 



dd 

+^cosC (- cosC + sinCcos (0 - (pp)) 



+ — sm(0 - (pp) 



6{ri-y) drid<j)dC. (40) 

This form is used in Eq. (I33]l to calculate the transverse-momentum spectra. 
The interesting feature of this formula is that the derivative of the distance d 
with respect to the spacetime rapidity rj has disappeared. This means that 
for each value of the rapidity (and spacetime rapidity) one calculates the 
spectra using the function d defined exactly for this value of rapidity. Such 
a function is provided by the hydrodynamic code which should be executed 
also for exactly the same value of rapidity (treated as a parameter). We 
thus see that the Cooper-Frye prescription is consistent with the physical 
picture of non-interacting transversally expanding 2D clusters. Moreover, 
the system under study is not necessarily boost-invariant. 



7. Results 

In Fig. |4]we show the PHENIX pion spectra for the centrality classes 20 
- 30% and 60 - 70% p5|, while in Fig. [5] we show the PHENIX data on the 
elliptic flow for the centrality classes 20 - 40% and 40 - 60% [26]. In both 
cases the data (solid lines) are compared with our hydrodynamic calculations 
(dashed lines) done for the appropriate values of the impact parameter. 
Following the PHENIX study performed within the Glauber model |25] we 
used the values: b = 7.25 fm for c = 20 - 30%, b = 7.94 fm for c = 20 - 40%, 
b = 10.26 fm for c = 40 - 60%, and b = 11.69 fm for c = 60 - 70%. The 
experimental spectrum of pions shown in Fig. [4] is the spectrum of positive 
pions multiplied by a factor of 3. This is done to account for the total hadron 
multiplicity. 

Our results were obtained for the initial central temperature Tj = 250 
MeV and for four different values of the final temperature Tf = 200, 180, 
160 and 140 MeV. Similarly to the results obtained by Heinz and Wong we 
observe that a lower than usual (i.e., lower than for 3D expansion) initial 
temperature is required to describe correctly the slope of the experimental 
spectrum. This effect is related to the generation of stronger transverse flow 
in the case of 2D expansion. In agreement with Ref. [20] we also find that 
the shape of the spectrum is quite insensitive to the final temperature since 
the lowering of the temperature in the distribution function is compensated 
by the increase of the magnitude of the transverse flow. This property can be 
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understood in our case as an effect of the exact conservation of the transverse 
energy and entropy. Since the entropy current is proportional to the particle 
current, the particle number is also conserved. Combing the conservation 
laws with the fact that particles have fixed rapidity we find that the average 
transverse momentum must be exactly conserved 

— dE 

-^ = ^ = (P^) = const. (41) 
dy 

The conservations of the average transverse momentum and the particle 
number do not allow for substantial changes of the slope of the spectrum. 
This behavior is shown in Fig. [4], where the spectra for different final tem- 
peratures are shown to be very much similar. 

The results presented in Fig. |4] were obtained with the normalization 
factor no = 1. The correct normalization for the peripheral collisions may 
be obtained by the small reduction of reo, while the correct normalization for 
the central collisions requires larger values of no. Dependence of the normal- 
ization factor on the centrality reflects the experimental fact that the total 
multiplicity grows faster with centrality than the number of the wounded 
nucleons. Frequently, this effect is understood as the extra contribution from 
the binary collisions. A possible explanation of such a centrality dependence 
is also offered by the model of wounded quarks and diquarks [SllEH]- Since 
such effects are not included in the form of our initial condition ((28]) , different 
values of no should be applied for different values of the centrality. 

In Fig. Owe show our main results concerning the elliptic flow coefficient. 
The very striking observation is that the elliptic flow is large and for lower 
values of the final temperature it even exceeds the data. The origin of 
this effect is the genuine two-dimensional hydrodynamic expansion which 
transforms the initial thermal energy exclusively into the transverse flow. We 
observe that for peripheral collisions the elliptic flow coefficient V2 becomes 
compatible with the data already at the temperature as high as 200 MeV. 
For more central collisions the experimental values are recovered at smaller 
but still high values of Tf . 



8. Conclusions 

In this paper we have shown that the idea that the system created in 
ultra-relativistic heavy-ion collisions undergoes the hydrodynamic expan- 
sion only in the transverse direction is compatible with the experimental 
data. With a suitable choice of the initial and final temperatures one is 
able to describe the measured hadron spectra and the elliptic flow coeffi- 
cient V2. Clearly, further investigations and developments of such a simple 
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idea should be performed to address more subtle aspects of liadron produc- 
tion. In particular, the present description should be supplemented with the 
hadronization model. 

Acknowledgments: We thank Professor Andrzej Bialas for his inspirations 
to continue this research and for the critical reading of the manuscript. 



9. Appendix 

The standard definition of the energy-momentum tensor as the second 
moment of the distribution function in the momentum space yields 

2ttt 

where (sinhx = sh x, coshx = ch x) 

/ da^T] (2uq + u^) chrjSuoUx chrjSuoUy chr/shr/ (2uq + n^) \ 



chr]3uoUx {1 + 3uxUx) 'iUxUy shr]3uoUx 

ch.r]3uoUy 3uyUx {1 + 3uyUy) sh.r]3uoUy 

\ chr/shry (2'Ug + shr]3uoUx shr]3uQUy sh^ rj (2uq + u\) ) 



This expression may be rewritten in the compact form, presented in the 
lower line in Eq. ((HI), if the definitions of the four- vectors and are 
used. 
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Fig. 4. Transverse-momentum spectra for the centrality c = 20-30% (upper part) 
and c = 60-70% (lower part). The PHENIX experimental results [25] (solid lines) 
are compared to the model calculations (dashed lines). The initial central tem- 
perature, Ti = 250 MeV, is the same in all cases. The four different values of the 
final temperature are considered: Tf = 200, 180, 160 and 140 MeV. As discussed 
in the text in more detail, the spectra are to large extent independent of the fi- 
nal temperature (the dashed fines overlap). The normalization factor no = 1 was 
used. 
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Fig. 5. The elliptic flow coefficient V2 for the centrality c = 20-40% (upper part) 
and c = 40-60% (lower part). The soHd lines represent the PHENIX data for pions 
and kaons [26] (the width indicates the experimental error) . The dashed lines show 
our model results. The values of the initial and final temperatures are the same as 
in Fig. 4. 



